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In this talk we report on a recent lattice investigation of the Landau gauge gluon and ghost propagators in 
pure SU{S) lattice gauge theory with a special emphasis on the Gribov copy problem. In the (infrared) region 
of momenta < 0.3 GeV'^ we find the corresponding MOM scheme running coupling Qs(q^) to rise in q. 
We also report on a first SU{3) computation of the ghost-gluon vertex function showing that it deviates only 
weakly from being constant. In addition we study the spectrum of low-lying eigenvalues and eigenfunctions of 
the Faddeev-Popov operator as well as the spectral representation of the ghost propagator. 



1. INTRODUCTION 

The infrared behavior of the Landau gauge 
gluon and ghost propagators in QCD is closely re- 
lated to the confinement phenomenon formulated 
in terms of the Gribov- Zwanziger horizon condi- 
tion and the Kugo-Ojima criterion Assum- 
ing a constant ghost-gluon vertex function the 
propagators also provide a nonperturbative deter- 
mination of the running QCD coupHng as(g^) 
in a MOM scheme. They are a subject of inten- 
sive studies within the framework of truncated 
systems of Dyson- Schwinger (DS) equations and 
within the lattice approach. 

In the infinite 4- volume limit the DS approach 
has led to an intertwined infrared power behavior 
of the gluon and ghost dressing functions |3I4I5| 



tor G{q'^). As a consequence, the MOM scheme 
running coupling 



in 



(2) 



Zn{q') ^ q'D{q^)o,{q^f\ 
Zciq') ^ q^G{q^)^{qT^ 



(1) 



with the same k « 0.59 , i.e. a vanishing gluon 
propagator D(q^) in the infrared occurs in inti- 
mate connection with a diverging ghost propaga- 
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is turning to a finite fix point in the infrared limit 
As we shall show below, our lattice simula- 
tions do not agree with this expectation. as{q^) 
seems to tend to zero in this limit. 

Possible solutions of the discrepancy have been 
discussed from different points of view in the re- 
cent literature. In |Q the DS approach has been 
studied on a finite 4-torus with the same trun- 
cated set of equations as for the infinite volume. 
as(g^) was shown to tend to zero for — > in 
one-to-one correspondence with what one finds on 
the lattice. This would indicate very strong finite- 
size effects and a slow convergence to the infinite- 
volume limit. However, we do not find any other 
indication for such a strong finite-size effect in 
our lattice simulations, except the convergence to 
the infinite- volume limit would be extremely slow. 
An alternative resolution of the problem has been 
proposed by Boucaud et al. [7] . These authors ar- 
gue the ghost-gluon vertex in the infrared might 
contain g^-dependent contributions which could 
modify the DS results for the mentioned propa- 
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gators^. However, recent detailed DS studies of 
the ghost-gluon vertex did not provide hints for 
such a modification 8 . Thus, at present there 
is no solution of the puzzle. In any case more 
thorough lattice investigations are desirable. 

For SU{2) extensive lattice investigations can 
be found in (SJ. Unfortunately, the authors did 
not reach the interesting infrared zone where the 
mentioned inconsistencies become visible. The 
same holds for SU{'i) lattice computations of 
the ghost and giuon propagators as reported in 
We have been pursuing an analogous 
study for the S'C/(3) case with the special em- 
phasis of the Gribov copy problem More- 
over, we have investigated the spectral properties 
of the Faddeev-Popov operator The low- 

lying eigenmodes of the latter are expected to 
be intimately related to a diverging ghost prop- 
agator. In addition, in ^1] we have reported on 
a first 5'[/(3) lattice computation of the ghost- 
ghost-gluon vertex at zero gluon momentum. We 
confirm, what has been found already for SU{2), 
namely the data are quite consistent with a con- 
stant vertex function 

In this talk we will give a short review of these 
investigations. In the meantime larger lattice 
sizes and better statistics have been reached con- 
firming our previous computations and conclu- 
sions. 

2. GHOST AND GLUON DRESSING 
FUNCTIONS 

We have simultaneously studied the gluon 
and ghost propagators in the quenched ap- 
proximation. SU{3) gauge field configurations 
U — {Ux,fj,} thermalized with the standard Wil- 
son gauge action have been put into the Landau 
gauge by iteratively maximizing the gauge func- 
tional 

1 4 

with gx G SU{S). It has numerous local maxima 
(Gribov copies), each satisfying the lattice Lan- 
dau gauge condition = 0. To explore to 
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what extent this ambiguity has a significant in- 
fluence on gauge dependent observables, we have 
gauge-fixed each thermalized configuration sev- 
eral times using the over-relaxation algorithm, 
starting from random gauge copies. For each con- 
figuration U, we have selected the first (f c) and 
the best (be) gauge copy (that with the largest 
functional value) for subsequent measurements. 
For details we refer to pTSj. 

It turns out that the Gribov ambiguity has no 
systematic influence on the infrared behavior of 
the gluon propagator. This holds as long as one 
restricts to periodic gauge transformations on the 
four-torus. In 16 the gauge orbits have been ex- 
tended to the full gauge symmetry including also 
non-periodic Z(iV) transformations. As a con- 
sequence the Gribov problem appears to be en- 
hanced even for the gluon propagator. Here we 
restrict ourselves to the standard case of periodic 
gauge transformations. In marked contrast to the 
weak Gribov copy dependence of the gluon prop- 
agator the ghost propagator at lower momenta 
depends on the selection of gauge copies. (For 
a study of the influence of Gribov copies on the 
ghost propagator in the SU (2) case see QT].) The 
effect of the Gribov copies is illustrated in Fig. ^ 
There we have plotted the fc - to - be ratios of the 
ghost and gluon dressing functions. Obviously 
for the ghost propagator the Gribov problem can 
cause 0(5%) deviations in the low-momentum re- 
gion {q < 1 GeV). However, a closer inspection of 
the data for the ghost propagator indicates that 
the influence of Gribov copies becomes weaker 
for increasing lattice size. This is in agreement 
with a recent claim by Zwanziger according to 
which in the infinite volume limit averaging over 
gauge copies in the Gribov region should lead to 
the same result as over copies restricted to the 
fundamental modular region JH|- In QHI simi- 
lar indications have been found for SU{2). Fig. El 
shows the ghost and gluon dressing functions ver- 
sus q"^ . In contrast to our previous papers |12I14| 
we have plotted only f e-copy results, where we 
have included some new data for lattices 24^ x 48 
and 32'* at /3 = 6.0. Both propagators have been 
renormalized separately for each /3 with the nor- 
malization condition Z— leLtq = fi^3 GeV. 
Note the deviation of the data for the asymmet- 
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Figure 1. The ratios Z^'^/Z^'^ for the dressing 
functions Zq (upper panel) and Zd (lower 
panel) determined on first (fc) and best (be) 
gauge copies. 



ric lattice 2A^ x 48. In principle, for better gauge 
copies the curve can become only less singular. 
Still it is very difficult to extract a possible power 
behavior at low in comparison with Eq. Q . 

3. THE GHOST-GLUON VERTEX AND 
THE RUNNING COUPLING 

In Fig. 13 we present the combined result for 
the running coupling Fits to the 1-loop and 
2-loop running coupling are also shown. Our 
new data points confirm our previous observa- 
tion that the running coupling monotonously de- 
creases with decreasing momentum in the range 

< 0.4 GeV^. We have plotted f c copy results, 
showing that this behavior is not a matter of se- 
lecting better gauge copies. The nice coincidence 
of the data for various lattice sizes (except the 
points refering to the asymmetric lattice) makes 
it hard to see how finite-volume effects can be 
blamed for the disagreement with the DS contin- 
uum result Eq. J^l- Thus, the arguments as given 
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Figure 2. The ghost dressing function Zq 
(above) and the gluon dressing function Z^ (be- 
low) vs. momentum both measured on fc 
gauge copies. 



in '7' should be checked in order hopefully to re- 
solve the conflict. 

One can ask, whether the ghost-ghost-gluon 
vertex renormalization function Zi[q^) is really 
constant at lower momenta. A recent inves- 
tigation of this function defined at vanishing 
gluon momentum for the SU (2) case 1151 supports 
that Zi{q^) « 1 at least for momenta larger than 
1 GeV. We have performed an analogous study for 
Zi{q^) in the case of SU{?>) gluodynamics. Our 
(still preliminary) results are presented in Fig. 2| 
There is a slight variation visible in the interval 
0.3 GeV^ < < 5 GeV^. But this weak devia- 
tion from being constant will not have a dramatic 
influence on the running coupling. 
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Figure 3. The momentum dependence of the 
running coupUng asiq^) measured on fc gauge 
copies. 



4. LOW-LYING MODES OF THE 
FADDEEV-POPOV OPERATOR 

Now let us discuss the spectrum of the low- 
lying eigenvalues of the Faddeev-Popov (F-P) 
operator. For each maximum of the gauge func- 
tional Fu [g] , besides of its discarded eight trivial 
zero modes, the eigenvalues of the F-P operator 
are positive. The corresponding gauge field con- 
figuration is said to lie within the Gribov region. 
If the lowest non-trivial eigenvalue tends to zero 
the configuration approaches the so-called Gribov 
horizon. 

Having the lowest n eigenvalues and eigenvec- 
tors determined, one can construct an approxi- 
mation to the ghost propagator 

n 

with ^i{k) being the Fourier transform of the 
«-th eigenmode (k denotes the lattice momen- 
tum). If all n = 8V eigenvalues and eigenvectors 
were known, the ghost propagator G{q^) would 
be completely determined. In the recent liter- 
ature jlHI2()j it is stated that the F-P operator 
acquires very small eigenvalues in the presence of 
vortex excitations. In any case, according to a 
popular belief it is an enhanced density of eigen- 



Figure 4. The inverse ghost-ghost-gluon vertex 
renormalization function Z^^{q'^) measured on f c 
gauge copies and normalized to unity at q = 3 
GeV. 



values near zero which causes the ghost propaga- 
tor to diverge stronger than l/g^ near zero mo- 
mentum |18j . 

At /3 = 6.2 we have extracted the 200 lowest 
(non-trival) eigenvalues and their corresponding 
eigenfunctions for lattice sizes 12'*, 16* as well as 
50 modes for 24^. Additionally, 90 eigenvalues 
have been extracted on a 24* lattice at = 5.8. 
This allows us to check how the low-lying eigen- 
values are shifted towards A = as the physical 
volume is increased. In order to clarify how the 
choice of gauge copy influences the spectrum, the 
eigenvalues and eigenvectors have been extracted 
separately on our sets of fc and be gauge-fixed 
configurations. 

In Fig. the distributions of the lowest Ai and 
second lowest A2 eigenvalues of the F-P opera- 
tor are shown for different (physical) volumes. 
There /i(A, A -I- A) represents the average num- 
ber (per configuration) of eigenvalues found in 
the intervall [A, A + A]. Open bars refer to the 
distribution on f c gauge copies while full bars to 
that on be copies. From this figure it is quite 
obvious that both eigenvalues, Ai and A2, are 
shifted to lower values as the physical volume 
is enlarged. In conjunction the spread of A val- 
ues decreases. We have found that the center of 
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Figure 5. The distribution h{X) of the lowest (left 
panels) and second lowest (right panels) Faddeev- 
Popov eigenvalues Ai and A2. Filled columns refer 
to be gauge copies, while open columns refer to 
f c copies. 



those distributions tends towards zero faster than 
Here L refers to the linear lattice extension 
in physical units. For example, we have found 
(Ai)(L) oc l/L^+e with e « 0.16(4). It is also vis- 
ible that the eigenvalues Ai and A2 on fc gauge 
copies are on average lower than those obtained 
on be copies. 

Eq. Q suggests that low-lying eigenvalues and 
eigenvectors have a major impact on the ghost 
propagator at the lowest momenta. Actually, this 
is not so easy to predict. We have studied, to 
what extent the lowest n modes saturate the esti- 
mator of the ghost propagator G„(g^) on a given 
set of gauge copies. We show the result in Fig. 
for the lowest {qi ) and the second lowest momen- 
tum (52) available on lattices of different size at 
/3 — 6.2. The estimates for Gn{q^) have been di- 
vided by the full (conjugate gradient) estimator 
G{q^) in order to compare the saturation for dif- 
ferent volumes. Considering first the lowest mo- 
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Figure 6. The ratio of the ghost propagator 
Gn{(f) approximated in accordance with Eq. |@J 
to the corresponding full estimator G{q^) as a 
function of n for the lowest qi and second low- 
est momentum (72 {(3 = 6.2, lattice sizes 12^, 16'' 
and 24^). 



mentum qi we observe from Fig. |H1 that the rates 
of convergence differ, albeit slightly, for the three 
different lattice sizes. The relative deficit rises 
with the lattice volume. In order to reproduce the 
ghost propagator within a few percent, 150 .. . 200 
eigenmodes turn out to be sufficient. For the sec- 
ond lowest momentum 52 even 200 eigenmodes 
are far from saturating the result. 

5. CONCLUSIONS 

We have studied the low-momentum region of 
Landau gauge gluodynamics using Monte Carlo 
simulations with the Wilson plaquctte action and 
various lattice sizes from 16'' to 32^. For the in- 
verse bare coupling we haven chosen (3 ~ 5.8, 
6.0 and 6.2. In this way we have reached mo- 
menta down to q^ ~ O.lGeV^. We have presented 
data mainly for the f c gauge copies, but we have 
shown the ghost propagator to become less singu- 
lar within a 0(5%) deviation, when better gauge 
copies are taken. We have found indications that 
the Gribov effect weakens as the volume is in- 
creasing. Towards g ^ in the infrared mo- 



6 



mentum region, the gluon dressing function was 
shown to decrease, while the ghost dressing func- 
tion turned out to rise. However, the connected 
power laws predicted by the infinite-volume DS 
approach cannot be confirmed from our data. 
Correspondingly, the behavior of the running cou- 
pling as{q^) in a suitable momentum subtraction 
scheme (based on the ghost-ghost-gluon vertex) 
does not approach the expected finite limit. In- 
stead, this coupling has been found to decrease 
for lower momenta after passing a turnover at 
g2 «0.4 GcV2. 

In addition, we have investigated the spectral 
properties of the F-P operator and its relation 
to the ghost propagator for various lattice sizes. 
As argued in |18| we have found the low-lying 
eigenvalues to be shifted towards zero (the con- 
figurations approaching the Gribov horizon) as 
the volume is increased. The low-lying eigen- 
values extracted on be gauge copies (those with 
the largest functional value) are larger on average 
than those on f c copies. Thus, for finite volumes 
better gauge-fixing (in terms of the gauge func- 
tional) keeps the configurations at larger distance 
from the Gribov horizon. The study of the contri- 
butions to the ghost propagator coming from the 
lowest eigenmodes shows that even at the lowest 
momentum one needs a large number of modes in 
order to saturate the propagator. For increasing 
volumes and larger lattice momenta the conver- 
gence of the eigenmode expansion becomes even 
slower. 

All simulations have been done on the IBM pSeries 
690 at HLRN. We are grateful to Hinnerk Stiiben for 
providing parts of his code. One of us (M. M.-P.) ac- 
knowledges useful discussions with V.K. Mitrjushkin. 
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